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Dilution effects on spin-1/2 quantum Heisenberg antiferromagnets with a non-magnetic 
spin-gapped ground state are studied by means of the qunatum Monte Carlo simulation. 

■ In the site-diluted system, an antiferromagnetic long-range order (AF-LRO) is induced at 
an infinitesimal concentration of dilution due to an effective coupling J mn between induced 

. magnetic moments. In the bond-diluted case, on the other hand, the AF-LRO is not induced 

up to a certain concentration of dilution due to singlet pairs reformed by an AF coupling J a f 

S between the edge spins of the diluted bond through the two-dimensional shortest paths. The 
competition between J mn and J a f yields peculiar phenomena in the bond-diluted system. 

: 

§1. Introduction 

o : 

l— ', Dilution effects on low-dimensional antiferromagnets with a spin-gapped ground 

' state exhibit various interesting phenomena inherent to low-dimensional quantum 

■ spin systems. In the first inorganic spin-Peierls compound CuGe03, the dimerized 
state is realized at low temperature and the system has the spin-gapped ground 

| state. x ) When a small amount of the non-magnetic impurity Zn or Mg is substituted 

for Cu, the spin-gapped ground state of the pure system changes to the antifer- 
romagnetic (AF) long-range ordered state, 2 ) " 4 ) thereby the lattice dimerization is 
preserved. The impurity-induced AF long-range order (LRO) has been observed 
also for other spin-gapped materials such as the two-leg ladder SrCu203 5 ** and the 
Haldane compound PbNi2V20s. 6 ) The non-magnetic-impurity-induced AF-LRO has 
been numerically investigated by simulating the site-diluted AF Heisenberg model 

*^ ■ with the spin-gapped ground state. 7 > ' 9 ' In the site-diluted system, the AF-LRO is 

induced even at 1 % of site dilution for the AF Heisenberg model of the coupled 
alternating chains, 9 ) suggesting the critical concentration of non-magnetic impuri- 
ties for AF-LRO to come up is actually (see also Fig. 1 below). When spins are 

. ,— i | removed from sites, the effective spins are induced, whose magnitudes are peaked 

at the neighboring site connected to the removed site by a strong bond. They have 
the same extent as the correlation length of the pure system. Furthermore, there 



exists an effective coupling J mn oc (— l)l T ' m rn+1 lexp[— l/^J £p£p] between the effec- 
tive spins, 10 )' n ) which preserves the staggerdness with respect to the original lattice 
and decays exponentially as the distance between the effective spins, I = \r m — r n \, 
increases. Here ^p' y are the anisotropic correlation lengths of the pure system in the 
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two directions on a square lattice. Thus the AF-LRO is induced at an infinitesimal 
concentration of dilution for the two-dimensional (2D) site-diluted system. 

In the present work we have investigated ground-state phase transitions in the 
bond-diluted 2D spin-1/2 AF Heisenberg model of the coupled alternating chains on 
a square lattice by means of the quantum Monte Carlo (QMC) simulation with the 
continuous-imaginary-time loop algorithm. In the bond-diluted case, in contrast to 
the site-diluted one, a finite critical concentration of dilution, x c , exists, above which 
the system has the AF-LRO. In the extremely- low-concentration region, the system 
has a spin gap independent of the concentration of the diluted bond, x. The spin gap 
A is proportional to the squared interchain coupling J' 2 significantly smaller than 
the original spin gap in the pure system, A p . This spin gap A is caused by singlet 
pairs of two edge spins on the diluted bond by the AF coupling J a f through the 2D 
shortest paths. Furthermore, when x increases, A starts to decrease, and tends to 
vanish at x which is smaller than x c . This may suggest the existence of the quantum 
Griffiths phase in this system. 



§2. Model and method 



The bond-diluted AF Heisenberg model of the coupled alternating chains on a 
square lattice investigated in the present work is described by the Hamiltonian 



H - e (2i,j)(2i+l,i) S 2ij • S 2i+l,j + « 51 e (2i+lj)(2i+2j)S2i+l,j ' S 2 i+2, 



h3 



hJ 



(2-1) 



where 1 and a (0 < a < 1) are the AF intrachain alternating coupling constants, 
J' (> 0) the AF interchain coupling constant, and Sjj the quantum spin operator 




Fig. 1. Concentration dependences of the staggered magnetization in the site- 9 - 1 and bond-diluted 
systems with a — J' — 0.5. In the bond-diluted case there is the finite critical concentration 
x c ~ 0.05. The percolation thresholds on the site and bond processes are denoted by Xp and 
Xp, respectively. All the lines are guides to eyes. 
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with S = 1/2 at site (i, j). Randomly quenched bond occupation factors {c<i,j)(k,l)} 
independently take either 1 or with probability 1 — x and x, respectively. The pure 
system described by (2-1) with £uj)(k,l) = 1 f° r an bonds is in either a spin-gapped 
phase or an AF long-range ordered phase at zero temperature depending on the 
strengths of a and J'. 12 ) In the a = 0.5 system which we examine in the present 
work, there exists the spin-gapped phase for J' < J' c ~ 0.55. 

The QMC simulations with the continuous-imaginary-time loop algorithm 13 ) " 15 ) 
are carried out on L x L square lattices with the periodic boundary conditions. For 
each sample with a bond-diluted configuration, 10 3 ~ 10 4 Monte Carlo steps (MCS) 
are spent for measurement after 500 ~ 10 3 MCS for thermalization. Sample average 
for dilution is taken over 10 ~ 10 3 samples depending on L, x, and the temperature 
T. 

§3. Results and discussions 

The staggered magnetization M s (x) at T = is evaluated by 

M 2 (x) = ^ Um , (34) 

where the bracket (• • •) denotes both the thermal and random averages. The concen- 
tration dependences of M s in the site- 9 ) and bond-diluted systems with a = J' = 0.5 
are shown in Fig. 1. In the site-diluted case, the value of M s becomes finite even at 
1 % of site dilution. This is caused by the strong correlation between the effective 
spins due to the effective coupling J mn mentioned in § 1. In the bond-diluted case, 
on the other hand, there is the finite critical concentration of x c ~ 0.05 even for 
a = J' = 0.5, which is the system near the critical point J' c ~ 0.55 for a = 0.5 and 
x = 0. Further dilution would yield the disappearance of M s near the percolation 
thresholds. 16 )' 17 ) 

Next we examine the spin gap A between the ground state and the first excited 
state to reveal the nature of the disordered phase below x c in the bond-diluted 
system. We adopt the second-moment method, by which the value of A^ 1 defined 
as 

(Z\( 2 ))- 1 = lim lim A / JEM ~, (3-2) 

is generally considered to be a good estimate of A. In Eq. (3-2) (3 = 1/T and 
S(ir, 7T, u) is the dynamic structure factor at momentum (n, n) described by 5(71", n, u) = 
Jq C(r)e* aJr dr, where C(r) is the correlation function of the staggered magnetization 
in the imaginary-time direction: 

W = 4gE /V(-1)I-^(t')S|(t' + t)) (3-3) 

with Sf{r) being the z component of the spin operator at the imaginary time r. Its 
x-dependence for the system with a = 0.5 and J' = 0.3 is shown in Fig. 3 below. 
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Fig. 2. Dependences of A (2) on x (= 1/2L 2 ) for a = 0.5 and L = 16 ~ 80. Points denote the QMC 
data for J' = 0.3, 0.25, 0.2, and 0.15 from top. The lines show the fitting function A (2) (x). In 
the inset the J' dependence of the spin gap is shown. The line with the slope of 1.5 is the guide 
to eyes. 



For a small x ( < 0.01) A^ rapidly decreases with increasing x. However as we will 
argue below this is an artifact of the second-moment method except for exactly at 
x = 0. 

Let us assume that, in the low-concentration limit {x — > 0), there exists a mid- 
gap A induced by dilution inside the original spin gap, A p , of the pure system. The 
former excitation mode contributes to C(r) in proportion to x so that it may less 
contribute to A^ as compared to the back ground mode of A p . More explicitly, 
C(t) in this situation and in the limit L — > oo may be written as 15 ) 



C(r) 



(3 A (3 (3A [3 

(1 — ax)(sinh — ^-^) _1 cos1i[(t — ^)^p] + ax(sinh-^-)~ 1 cosh[(r — -^)A\ , 

(3-4) 

where a is a non-trivial constant depending on a and J'. Substituting Eq. (3-4) for 
Eq. (3-2), we obtain the following relation between A^ 2 \ A p and A: 



AW{x) = A 



\ 



1 + 



l-ax A 



1 + 



I— ax / A 
ax v A D 



(3-5) 



To confirm that the above assumption is in fact the case, we simulate the systems in 
which only one strong bond is removed, i.e., x = 1/2L 2 . In Fig. 2 the x dependences 
of the QMC data of A^ are shown for J' = 0.3, 0.25, 0.2, 0.15 and a = 0.5. 
They are well fitted to Eq. (3-5), the lines in the figure. This implies that the 
spin gap A, or the mid- gap, is independent of x at least in the x-range shown in 
Fig. 2. By this fitting the value of A is obtained as shown in the inset of Fig. 2, 
in particular A = 0.138(5) for J' = 0.3 as indicated in Fig. 3. The existence of 
the spin gap proportional to J' 2 is interpreted that two edge spins induced at both 
ends of a diluted bond reform a singlet pair due to the effective coupling J a f (~ J' 2 ) 



a = 0.5, J' = 0.3 



1/f 
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Fig. 3. The concentration dependences of A^ corresponding to the upper bound of the spin gap 
and the inverse spatial correlation length for a — 0.5 and J' — 0.3. This system has the 
critical concentration x c = 0.11(1). The value of A^ at x — equals 0.31415(2). The dotted 
line describes the spin gap independent of x in the extremely-low-concentration region and the 
expected reduction in the low-concentration region. 



through the 2D shortest paths. 18 ) Since two edge spins on the diluted bond are on 
the different sublattices, the interaction is always AF. 

As x increases, the distinction between a mid-gap and the gap A p becomes 
obscure. Instead, the value of A^ gives us an upper bound of the spin gap A of 
our interest, whose x-dependence is schematically shown by the dotted line in Fig. 
3. Presumably A^ is a good estimate for A in the range x > 0.02. 

In Fig. 3 we also show the x dependence of the spatial correlation length £ x (ar) = 
lim^oo £ X (L; x) at T — > 0. The critical concentration x c is about 0.11, which is 
estimated by the finite-size-scaling analysis using the QMC data of £ x (L;x). Its 
details will be reported elsewhere. In marked contrast to l/£ x , A^ (~ A in the 
range x > 0.02) seems to vanish at a certain value x' c smaller than x c , which strongly 
suggests the existence of the quantum Griffiths phase at x' c < x < x c . Thus the 
bond-diluted system has the more peculiar properties than the site-diluted system. 

The x dependence of the Neel temperature in a bond-random system CuGei-^Si^Os 
has been investigated experimentally, 19 - ) ' 20 - ) and has turned out to be similar to that 
of the site-diluted system Cui_a;Mg x Ge03. Besides the problem whether the sub- 
stitution of Si for Ge simply corresponds to the bond dilution or not, we have to 
investigate the systems at finite but low temperatures, explicitly taking into account 
the interlayer coupling. Experimentally, analyses on systems with further smaller x 
and at lower temperatures are desirable. 



§4. Summary 



We investigated dilution effects on the spin-1/2 quantum AF Heisenberg model 
with the spin-gapped ground state. When the strong bonds are removed from the 
lattice for the spin-gapped ground state, two effective interactions with the couplings 
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J mn and J a f are induced between the effective spins. The competition between J mn 
and J a f leads to peculiar quantum phase transitions. In the low-concentration region, 
J a f between edge spins mainly contributes to magnetic properties and, therefore, the 
new spin gap due to singlet pairs reformed by J a f is induced within the original spin 
gap. The value of the spin gap decreases as x increases and tends to vanish at x' c 
which is less than x c or just at x c . If the former is the case, the quantum Griffiths 
phase would exist. Thus there are quite different magnetic properties at T = 
between the site- and bond-diluted systems. 
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